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We analyze theoretically polar molecules confined in planar arrays of one dimensional tubes. In 
the classical limit, if the number of tubes is finite, new types of "clustered Wigner crystals" with 
increasingly many molecules per unit cell can be stabilized by tuning the in-plane angle between the 
dipolar moments and the tube direction. Quantum mechanically, these phases melt into distinct 
"clustered Luttinger liquids." We calculate the phase diagram of the system and study the quantum 
melting of the clustered phases. We find that the requirements for exploring these phases are 
reachable in current experiments and discuss possible experimental signatures. 
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I. INTRODUCTION 

Systems with competing long-range interactions of- 
ten exhibit structures with emergent large length scales. 
Some examples include the formation of bubble and 
stripe domains in Langmuir-Blodgett films or in thin fer- 
romagnetic layers,!^ and the chain formation of magnetic 
particles in three-dimensional ferrofluids.'^- Long-range 
dipolar interactions in a back-gated t wo- dimensional elec- 
tron gas (2DEG) have been predicted^ to lead to the ex- 
istence of "microemulsion" phases intervening the Fermi 
liquid and the Wigner crystal phase. Similar microemul- 
sion phases may appear in 2DEGs subject to magnetic 
fields such that several Landau levels are occupied. -''S] 

Theoretically, quantum emulsion phases are challeng- 
ing to analyze since they involve structures at length 
scales ranging from the inter-particle distance to meso- 
scopic scales. In contrast, the conventional tools of many- 
body physics are mostly geared toward two particle cor- 
relations, such as paired states, magnetism and charge 
density wave. Experimentally, quantum emulsion phases 
are not easy to probe since transport measurements can 
only provide indirect evidence about their existence. Re- 
alizing long-range interactions with systems of cold polar 
moleculcs^'^^'^ can allow to explore emergent emulsion 
phases in a highly controllable setting. Moreover, in such 
systems few-body bound states,'^ trimer liquid phases,'^ 
and bound solitons^*' have been predicted. 

Here we demonstrate that the anisotropic and long- 
range character of dipolar interactions leads to new types 
of clustered crystal phases which appear at intermediate 
values of the interaction strength. Quantum mechani- 
cally, these phases melt into distinct "clustered Luttinger 
liquids" characterized by the decay of their density- 
density correlation functions. We calculate the phase di- 
agram and study the quantum melting of the clustered 
phases when tuning the orientation of the dipoles. Our 
calculations indicate that the clustered phases can be ex- 
plored under current experimental conditions. 
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FIG. 1. (Color online) Proposed setup to observe cluster for- 
mation of polar molecules (a). The classical phase diagram 
for two tubes as a function of the tilting angle and the par- 
ticle density p is shown in (b). Lobe shaped phases consisting 
of clusters with q particles per tube emerge. The phase sep- 
arated regions are indicated by the shaded layer surrounded 
with speckles. Inset: Optimized cluster configuration in one 
unit cell for — 1.01, upper graph, and (j) = 0.99, lower graph, 
at p = 0.7, corresponding to g = 2 and q — 3, respectively. 



II. EXPERIMENTAL SETUP 

We consider a setup in which polar molecules are con- 
fined to riT one-dimensional parallel tubes [see Fig.[l](a)], 
which can be realized by deep optical lattices.'^ The 
dipolar moments are aligned in the plane of the tubes 
at an angle (j) with respect to the tube direction. The 
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FIG. 2. (Color online) Phase diagram in the tilting angle (j) and the particle density p plane for (a) ut = 3, (b) ut ~ 4, and 
(c) Ut ~ 8 tubes. The phase separation regions are indicated by shaded layers. 



inter-tube distance 5 is used as unit of length throughout 
this work. 

The interaction energy between two molecules with 
dipolar moment m is 

where re^ is the inter-molecule displacement and m = 
/i(cos(/), sin(/), 0)"^. For tilting angles below the crit- 
ical angle (pc — arccos 1 / the interaction between 
molecules in the same tube is attractive and the sys- 
tem is unstable. Thus, we focus on dipolar orientations 
7r/2 > (j) > (pc where the intra-tube interaction is re- 
pulsive. Yet, molecules in different tubes attract when 
their displacement along the tubes is not too large. It is 
precisely this interplay between attraction and repulsion 
which leads to the formation of clusters. 



III. CLASSICAL LIMIT 

We first discuss the emergence of mesoscopic structures 
in the classical limit [Ti ^ Q). 

When the dipoles are oriented perpendicular to the 
tubes {(j) = 7r/2) the ground state is a Wigner crystal with 
TLT molecules per unit cell and periodicity where p 
is the linear density of molecules. Upon tilting the direc- 
tion of the dipoles toward the critical angle 4'c, phases 
with an increasingly complex unit cell are formed, be- 
fore eventually becoming unstable to collapse at (j) — (j)c- 
In these phases, the unit cell consist of q particles per 
tube forming a superlattice with periodicity qp^^. For 
example, a phase with ut = 2 and g = 3 is illustrated 
schematically in Fig.[l](a). 

We derived the phase diagram as a function of the 
tilting angle and the density p, by minimizing the clas- 
sical ground state energy with respect to the position of 
the molecules, allowing for arbitrary periodic structures 
with up to g = 8 molecules per unit cell in each tube. 
The phase diagram for = 2 is shown in Fig.Il|(b). 



At small densities we observe transitions to phases with 
monotonically increasing q when decreasing the tilt angle 
from 7r/2 toward (pc- Phases of a fixed value of q have 
a lobe-like structure, which bend with increasing density 
toward larger </>. Quite generally phases in Fig.[l](b) ter- 
minate by phase separated regions, indicated by a shaded 
layer surrounded with speckles. The phase separated re- 
gions are determined by the Maxwell construction which 
is applicable when the interfacial energy is positive.'^ It 
is possible that phases with g > 8, not captured by our 
present calculation, are favorable in some parts of the 
phase diagram. In particular, this is the case very close 
to (j) — (j)c where we find that q = 8 has the lowest energy. 

The origin of the cluster formation can be easily under- 
stood by considering the case (f> = pc- Then, the intra- 
tube repulsion is precisely zero. In order to maximize 
intcr-tube attraction, it is favorable to form a single clus- 
ter with a macroscopic number of particles, correspond- 
ing to a (J — ?> oo phase. As the angle is tuned toward (j)c^ 
there must be either an infinite sequence of transitions 
to increasingly higher values of g, or a macroscopically 
phase separated region.!^ 

Next, we discuss systems with more than two tubes. 
Results for riT = 3, 4, and 8 tubes are shown in Fig. [2] 
The phase diagrams for nx > 2 have a similar lobe struc- 
ture as in the ut = 2 case. The main difference is 
that with increasing number of tubes the lobes extend 
to higher values of the tilting angle (j). Thus, clustered 
phases might be easier to observe in systems with a larger 
number of tubes. As in the — 2 case, phase separated 
regions appear between phases of different q. 

A two-dimensional system which consists of an infinite 
number of tubes with dipoles aligned in the plane ex- 
hibits similar physics: the q = 1 Wigner crystal phase 
becomes locally unstable for (j) > pc- However, in this 
case, trial configurations with an increasing q have mono- 
tonically lower energy (we have tried structures with up 
to q = 128), indicating that the ground state may be 
phase separated. In the low density limit, the dipoles 
form infinitely long strings, which are mutually attrac- 
tive^" and thus one can show that the system is unstable 
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to macroscopic phase separation. Note that, for in-plane 
dipoles, the (logarithmicahy divergent) surface energy is 
posiiz?;epSI therefore macroscopic phase separation is pos- 
sible (unlike the out-of-plane case|31). 



IV. QUANTUM MECHANICAL ANALYSIS 

In quantum mechanical systems with continuous trans- 
lational symmetry, true long-range crystalline order ap- 
pears only in two dimensions or higher, even at zero 
temperature. In one-dimensional systems, the density- 
density correlations decay for large distances as a power 
law. Nevertheless, one can expect that upon melting the 
clustered Wigner crystal phases by quantum fluctuations, 
these phases will remain distinguishable by the nature of 
their quasi-long range correlations. We term the result- 
ing phases "clustered Luttinger liquids." 

In a clustered Luttinger liquid phase, the slowest- 
decaying component of the density-density correlations 
has a spatial period of A = qp~^ ■ In a bosonized descrip- 
tion, the fundamental harmonic of the density operator 
is therefore of the form cos[27r(a; -I- xo)p/q], where xq is a 
uniform shift of the crystalline configuration. In terms of 
the "counting field" defined relative to the crys- 

talline configuration we obtain for the bosonized density 
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where the dots represent higher harmonics. The factor 
q~^ present in the cosine alters the power law with which 
the density correlations decays: 
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Here, K is the Luttinger parameter. Microscopic 
consideration^^ suggest that K cc q. Therefore, the ex- 
ponent with which the density-density correlation func- 
tion decays is proportional to 1/q. Phases with larger q 
thus have a slower decay of the density-density correla- 
tion function, and are increasingly "classical" in nature. 

In order to make quantitative predictions about the 
phase diagram in the presence of quantum fluctuations, 
we have investigated a system of two tubes numeri- 
cally by means of density matrix renormalization group 
(DMRG)22123| simulations. To this end, we introduce the 
lattice Hamiltonian 



P^ V— \' 

TJ 2^ - j)a/6, a - /3]n«,jn/3j , (3) 



where ^ {c'^ j) destroys (creates) a particle at site i of 
tube a — 1,2, and n^^i = i^a i counts the number 




FIG. 3. (Color online) Particle density p{x) in the lower 
(red squares) and the upper (blue triangles) tube for dipolar 
strength 7 = 8, left column, and 7 = 50, right column, for 
a system of length L — 245, density p = 0.55~^ and lattice 
spacing a = 5/4. From top to bottom the tilt angle takes 
the values = {1.02, 0.99, 0.98, 0.97, 0.96} where for 7 = 
50 pronounced clusters with q = {1, 2, 3, 4, 6} are found, 
respectively. The data is obtained by DMRG for systems 
with open boundary conditions. 



of particles. Due to the strong on-site repulsion we treat 
the particles as hard-core, and therefore for the quantities 
we compute here (e.g., density distributions and ground- 
state energies), it does not matter whether the particles 
are bosons or fermions. The discrete Hamiltonian can 
represent a continuous system by taking the lattice spac- 
ing a — > while keeping the product ta^ = l/(2m) con- 
stant, where t is the hopping strength and m is the mass 
of the particles in the continuum. The primed sum in- 
dicates that the singular contribution where i — j and 
simultaneously a = (3 is omitted. The dipolar energy 
Vd{x,y) is given by Eq. ([ij with r = (x, y)^ . 

It is convenient to introduce the dimensionless quan- 
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FIG. 4. (Color online) Quantum phase diagram for a two 
tube system of length L — 245, density p — 0.55~^, and lattice 
spacing a = 5/4, as a function of the tilting angle (j) and 
the ratio between the kinetic and the interaction energy I/7. 
Right of the asterisks connected by lines cluster formation 
can be observed. To the left of this line, the ground state is 
a g = 1 Luttinger liquid. 



tity 7, which is the ratio between the typical dipolar 
interaction energy -Edip £^nd the typical kinetic energy 
-Ekin- These energies can be estimated as E^ip 



and Ekin ^ P respectively, and thus 7 ^ fi pm. In 
the limit of strong interactions, 7 3> 1, the system is ex- 
pected to be essentially classical, and the phase diagram 
is expected to be similar to that of Fig.jljwith the Wigner 
crystalline phases replaced by clustered Luttinger liquids. 
Conversely, for 7 <§; 1, quantum fluctuations dominate, 
and we expect only the q — 1 phase to survive. 

The particle density p{x) evaluated with DMRG for 
the Hamiltonian (3| exhibits clear signatures of clustered 
phases, see Fig. [3 We consider a system of finite length 
L — 24(5 with open boundary conditions and particle den- 
sity p = 0.56^^ and dipolar strength 7 = 50. The lattice 
constant is a = ^/4; no significant change in the results 
was observed when a was decreased to 6/6. Since the 
reflection symmetry about a plane perpendicular to the 
tube is broken for any tilting angle except </> = 7r/2, the 
density of the upper pu and the lower pi tube are slightly 
shifted. Additionally, we observe that the height of the 
peaks in the density decreases toward the trap center, 
consistent with quasi-long range order. Remarkably, the 
rate of the decay decreases strongly with increasing q, as 
expected from the Luttinger liquid analysis, Eq. ([2|. 

The complete quantum phase diagram for the two tube 
system as a function of the tilting angle and the ratio 
between kinetic and interaction energy I/7 is shown in 
Fig. [4] For 1/7 — 0, the results were obtained by classi- 
cal minimization of the interaction energy. The DMRG 
simulations are used to extend the results to I/7 > 0. 
The phases are determined from the density distribution 
by calculating the number of particles localized within 



one cluster. The clustered Luttinger liquid phases with 
q > 1 extend to considerably large values of I/7, making 
the realization of clustered phases feasible in experiments 
with cold dipolar molecules. 



EXPERIMENTAL IMPLICATIONS 



For typical densities of p = 10'*cm 



7 



{0.7, 3.3, 6.8, 49.4, 63.3} can be achieved in experiments 
with KRb, RbCs, NaK, NaCs, and LiCs, respectively. 
(See [201 for details on this calculation.) 

The density regime which is most favorable for ob- 
serving clustered phases is p ^ 0.5^~^ — 5~^. For 
p ^ lO^cm^^, this corresponds to an inter-tube separa- 
tion of 5 ~ 5 • 10~^cm, easily attainable using an optical 
potential created by a laser with wavelength ~ 1pm. 

NaCs and LiCs are thus the most promising candidates 
to realize clustered phases, due to their large dipolar mo- 
ments. In order to make the clustered phases more ro- 
bust, one can add a shallow periodic potential along the 
tubes. Such a periodic potential quenches the kinetic 
energy, thus increasing the effective value of 7. As a con- 
sequence cluster formation arises at much weaker dipolar 
moment also attainable by KRb. 

Other effects that can be important for experiments 
are (i) the incommensurability of the particle number 
with the cluster size, (ii) the shallow trap potential along 
the tube direction, (iii) the strong but finite transverse 
confinement, (iv) quantum fluctuations in the orienta- 
tion of the dipoles, and (v) finite temperature effects. 
Cluster formation is extremely stable with respect to 
(i) and (ii).'^ Incommensurability leads to a slight re- 
arrangement of clusters and the consequence of the shal- 
low trap along the tubes is merely that the distance be- 
tween the clusters is reduced. One-dimensional tubes 
are realized by strong transverse confinement potentials 
(iii). Therefore, we consider the interactions computed 
for molecules with transverse wavefunctions, correspond- 
ing to a parabolic confinement, and compare them to 
the bare, one-dimensional interactions. The renormal- 
ization of the interactions due to the transverse con- 
finement can be evaluated from a multipole expansion 
yielding AE^^^f E^ip ^ cr^/<5^, where a is the spread 
of the wave function in transverse direction-I^Sl Under 
standard experimental conditions a ^ 25nm and thus 
^-^dipZ-^dip ^ 0.0025. This ratio has to be compared 
with the relative energy difference between the clustered 
q > I and the uniform q = 1 phase, which typically is 
0.2 — 0.5 ^ AEj-jp/Edip. Therefore, the renormaliza- 
tion of the interaction energy due to the finite strength 
of the transverse confinement is no obstacle for the ob- 
servability of the clustered phases. The quantum fluc- 
tuations of the dipoles around the orientation of the 
electric field (iv) renormalize the dipolar potential by 
AE^ip/Edip < Re/S - 0.001, where Re is the bond 
length of the molecule. Thus, this effect is also small. 
The temperature scale (v) below which we expect strong 
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tendency toward cluster formation is proportional to the 
dipolar energy E^ip. In units of the Fermi temperature 
Tp ^ E]^ia the crossover temperature is Tcross — ctj, 
where the proportionality constant a can be estimated 
from the relative energy difference of the uniform and 
the clustered phases, i.e., a ^ 0.2 — 0.5. Depending on 
the molecule^Sl the crossover temperature is Tdoss ^ Tp . 

VI. CONCLUSIONS AND OUTLOOK 

In summary, the formation of clustered phases of po- 
lar molecules is a distinct consequence of the long-range 
and anisotropic nature of their interactions. We find 
that the clustered phases can be explored under cur- 
rent experimental conditions. A variety of techniques 
can be employed to observe the clustered phases, includ- 
ing elastic light scattering,'^^! noise correlations in time- 
of-flight i mage s,'^ and optical quantum nondemolition 
detection.'i^'^Sl Cluster formation should also contribute 
additional dissociation resonances in lattice modulation 
experiments^^ and rf spectroscopy.f^SlEni 

A particular exciting direction for future research is to 
study excitations in this system. For the simple q — 1 
phase, the unit cell consists of a single particle per tube 
and thus only the acoustic mode exists. However, in 
clustered crystal phases with q > 1 also optical branches 
should exist in the excitation spectra. 
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SUPPLEMENTARY MATERIAL 

Incommensurability and longitudinal confinement. — 
Here, we investigate in detail the effects of the particle 
number incommensurability with the cluster size and of 
the shallow confining potential along the tube direction. 

To demonstrate that the cluster formation is stable 
against the incommensurability of the particle number 
with the cluster size, we evaluate the density distribu- 
tion for systems with = 10, = 12, and A^ = 14 polar 
molecules per tube, respectively. The remaining system 
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FIG. 5. (Color online) Particle density p{x) in the lower (red 
squares) and in the upper (blue triangles) tube for the ratio 
of the interaction to kinetic energy 7 — 2.4, tilt angle 0=1, 
length L — 805, and lattice constant a — S for systems with 
(a) TV = 10, (b) N = 12, and (c) A = 14 particles per tube. 
The cluster formation is insensitive as to particle densities 
which are incommensurable with the cluster size. 



parameters are: ratio of interaction and the kinetic en- 
ergy 7 — 2.4, tilt angle (f> — 1, length L — 806, and lattice 
constant a = S. The cluster formation is commensurable 
with A = 12 particles, where four clusters with q = 3 
particles are found, see Fig.[5](b). However, both A = 10 
and A = 14 are incommensurable with the q = 3 phase. 
Regardless, we observe a pronounced cluster formation, 
see Fig.[5](a) and (c). 

Next, we study the effects of the harmonic confining 
potential along the tube direction, which we incorpo- 
rated into the Hamiltonian by a variable on-site energy 
ej = C(ja — L/2)^/L^. It turns out that the formation of 
clusters is also very insensitive to the confining potential. 
Results for the particle distribution are shown in Fig. [6] 
for (a) C = O.lt, (b) C = l.Ot, and (c) C = 4.0t. For the 
considered values of C, the main effect of the harmonic 
confinement potential is that the distance between clus- 
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FIG. 6. (Color online) Particle density p{x) for the same system as described in the caption of Fig. [5] but with A'' = 12 
throughout and a confining potential of strength (a) C = O.lt, (b) C = l.Qt, and (c) C — 4.0t. The cluster formation is robust 
with respect to the trap potential. The mere effect of the trap in the studied regime of the confining strength is that the 
distance between the clusters is reduced. 



ters of particles is reduced as compared to the uniform 
case shown in Fig.[5](b). 

Multipole expansion. — The renormalization of the in- 
teractions due to the fluctuations of the dipoles around 
the orientation of the applied electric field and due to 
the finite spread of the transverse wavcfunctions, can be 
evaluated from a multipole expansion of the electric po- 
tential 



form Ee2 = J^^eiSex) 



21 + 1 r'+i 



(4) 



where qi„i — J Yi*^r^ p{r)(Pr with the charge distribution 

The largest effect of this renormalization concerns 
dipoles which are separated roughly by S. The electric 
field induced from one dipole at the location of the second 
dipole is 

E(5e^) = V<j){Sex) , 

and the potential energy is given by = mE{5ex). The 
leading order term is clearly the dipole-dipole interaction 
E22 = mE2((5e:r)- 

The next order is the dipole-quadrupole interaction 
i?42, induced from the fluctuations of the dipoles around 
the electric field. This contribution is nonvanishing, 
since heteronuclear polar molecules are considered. From 
the multipole expansion Q we infer the microscopic 
quadrupolar renormalization 



E42 
E22 



const. 



Re 



with const. ^ 1 and Re the bond length of the molecule. 

As the transverse confinement obeys rotational sym- 
metry, there is no quadrupole contribution from the 
transverse wavefunction. Therefore, the first non- 
vanishing contribution due to the spread of the trans- 
verse wavefunction is of the dipole-seipole (6 charges) 



E22 ^ 



(5) 



where we used p+(r) = /5-(r) — q/ {na^) exp[—{y'^ + 
z^)/a-^]S{x) and thus qim ^ qReCr''^^. 

Local stability of clustered phases. — One concern could 
be that the uniform q = 1 crystalline phase may be 
metastable, and therefore the system might not be able to 
find its ground state on experimental timescales. Thus, 
it is useful to explicitly evaluate the stability of the uni- 
form q = 1 phase. A system is locally unstable if its 
compressibility k is negative. The compressibility can be 
evaluated from — d^e/d{p~^)^ , where e is the energy 
per particle and p the particle density. Thus it is useful 
to evaluate the spinodal line, which separates the regions 
with positive and negative compressibility, see asterisks 
connected by lines in Fig.[7](a). On the right hand side of 
the spinodal line, close to the critical angle (j)^, there is a 
large region where the uniform q — 1 phase is locally un- 
stable. In fact, it can be shown that any clustered phase 
(with an arbitrary q) must become locally unstable for a 
sufficiently large density and sufficiently close to (j) = (j)c- 

Assume that (j) = 'Pc-, such that the intra-tube inter- 
action is zero. In the high density limit, the energy per 
particle is proportional to the particle density times the 
integrated inter-tube energy. Explicitly integrating the 
inter-tube interaction gives e — — p? p siv? {(p) / 5"^ . There- 
fore kT^ < 0. As we move away from the critical angle, 
the repulsive intra-tube interaction gives a contribution 
proportional to {4> — 4>c)p^ ■ Adding these two contribu- 
tions and calculating the compressibility, we get that the 
critical density p*{4>*) at which the phase becomes un- 
stable satisfies p* cx (0* — (/'c)^- 

This argument captures the main characteristics of the 
spinodal line at high densities. Below a certain density, 
the spinodal line terminates and the phase becomes lo- 
cally stable, as can be seen in Fig.[7](a) for q = 1. In 
Fig.[7](b) we show the calculated ground state energy of 
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the 9 = 1 phase for (p = 0.99, showing the region in which 
K < 0. The relative extend and shape of the spinodal hne 
for systems with a larger but finite number of tubes, as 
shown in Fig.[2j is similar to the case of two tubes. 

At the critical angle (/> = i^c the ground state consists 
of dipoles which are arbitrarily close to each other cor- 
responding to a g = oo phase. From this and from the 
analysis of the compressibility follows that, there must 
be either an infinite sequence of phases with increasing 
q, or a region of macroscopic phase separation between a 
phase of finite q and vacuum upon approaching (j) = <f>c- 

Infinite number of tubes. — Next, we study the phase 
diagram in the case of an infinite number of tubes. We 
find that the q — 1 phase is locally unstable in a re- 
gion 7r/2 > (j) > (pc- Our results are consistent with 
macroscopic phase separation between the q = 1 phase 
and vacuum. So unlike the case with a finite number of 
tubes, intermediate phases with q > 1 are not realized. 

The classical ground state energy of the two- 
dimensional system is evaluated by Ewald summation 
techniques .'SiHlIl'yye allow for a unit cell containing q par- 
ticles. The unit cell is a parallelogram of length 
height (5, and an arbitrary angle. In addition, for q > 1, 
we assume that the particles within a unit cell are equally 
spaced. The angle and the intra-unit cell spacing are 
treated as variational parameters. The same scheme was 
used to optimize the ground state energy in the case of 
a finite number of tubes, and gave excellent results com- 
pared to a full optimization with respect to the positions 
of all the particles in the unit cell presented in the main 
text [see also inset of Fig.[T](b)]. 

Fig. [8] shows the spinodal line for the q = 1 phase, at 
which the compressibility is zero. To the right of this line, 
the compressibility is negative, and the q ~ 1 phase is 
locally unstable. In the unstable regime, phases with q > 
1 are lower in energy than the q — 1 phase; however, the 
energy decreases monotonically with increasing g, up to 
the largest q that we tried (g = 128), and is always higher 
than a phase separated state with macroscopic regions 
of g = 1 and vacuum. In the low density limit strings 
of polar molecules are formed, whose mutual interaction 
is attractive. Therefore, at sufficiently low densities the 
system always phase separates. 

In two-dimensional systems with dipolar interactions, 
the surface tension between two phases of different den- 
sity diverges logarithmically.Sl For dipoles pointing out 
of the plane, the divergent term is negative, and hence 
macroscopic phase separation is always unstable toward 
forming "microemulsion" phases with emergent, meso- 
scopic structures. For in-plane dipoles, however, the 
divergent energy is positive, and therefore macroscopic 
phase separation is possible. 

To demonstrate this, we consider a system with a linear 
domain wall between two phases A and B. The dipolar 
contribution to the total energy is 
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FIG. 7. (Color online) (a) Phase diagram of dipolar 
molecules confined in two tubes as in Fig. [T] and spinodal line 
for the g = 1 phase, asterisks connected by lines, which sep- 
arates regions of positive and negative compressibiHty. (b) 
Ground state energy of the g = 1 phase for (j) — 0.99. In the 
region between the dashed vertical lines the curvature of the 
energy per particle is negative indicating the density regime 
where the uniform phase is certainly unstable. 
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FIG. 8. (Color online) Spinodal line for a system with in- 
finitely many tubes, which indicates the region of instability 
of the q — 1 Wigner crystal phase. 
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where p{r) is the density at position r, which is pA in 
phase A and ps in phase B. Alternatively, the energy 
can be written as 

E ^ Vaca + Vbcb + , (7) 

where eA(e_B) is the energy density of the homogeneous 
A{B) phase, Va{Vb) is the area of phase A{B), and Ecr 
the surface energy between the two phases. The bulk 
contributions are defined as 

VAiB)eA{B) = \ f dr j dr' p\g^V{v-v') , 

where ^a{B) defines the domain of phase A{B). The inte- 
gral over r' is taken over the entire plane . Comparing 
Eqs. Q and ([t]) yields the surface energy 

Ea^-\ I dr f dv'{pA - PB?V{v - r') . 

For dipolar interactions, E^/L (where L is the linear di- 
mension of the system) diverges as log(L) as L — > oo. In 
the following, we determine the sign of E^ for a straight 
domain wall oriented along the y-direction. In order to 
make contact to the results of Ref. [H where the dipoles 
are pointing out of plane, we consider a general dipole 
orientation: 

m = p{cos(j)cos6, sin cos 61, sin 6*)"^ , 

where the angle 9 is measured from the plane. An explicit 
calculation reveals that the surface energy vanishes at an 
angle 9c given by 

9c = arccos ^ / — . 

V 1 -I- COS"' 

For 9 > 9c, Ecr < 0, and macroscopic phase separation is 
impossible; ior 9 < Oc, Ecr > Q and phase separation is 
allowed. In particular, for 9 — 7r/2, Ec, < (consistently 
with Ref. 4), whereas for 9 — 0, Ecr > for any value of 

0. ^ 

In summary, in the case of a two-dimensional array 
of tubes, it seems that the q = 1 phase terminates at a 
first-order transition to vacuum, and no higher q phases 
exist. 

Numerical values for 7. — The numerical data used to 
evaluate the ratio between the dipolar and the kinetic 
energy 7 is summarized in Tab.|l] 

Luttinger parameter K . — In order to reveal the depen- 
dence of the Luttinger parameter K of the clustered sys- 
tem on q, we relate it to the Luttinger parameter Kq of 
the uniform {q — 1) system. The low-energy properties 
of a 5 = 1 liquid of polar molecules in one-dimensional 



tubes are described bjEHSH 

TABLE I. Numerical values for the ratio between the dipo- 
lar and the kinetic energy 7 for dipoles AB and the quan- 
tities which are necessary to estimate this ratio. The dipole 
strength is denoted as p and mA {niB) is the mass of con- 
stituent A (B). We assume the linear density to be p = 
10*cm-\ 



AB 


p/De 


mA/u 


mB /u 


7 


KRlP 


0.6 


39.1 


85.5 


0.7 


RbCiP 


1.0 


85.5 


132.9 


3.3 


NaKP 


2.7 


23.0 


39.1 


6.8 


NaCsP 


4.6 


23.0 


132.9 


49.4 


LiCpl 


5.5 


6.9 


132.9 


63.3 



where we introduced the microscopic Luttinger parame- 
ters Kq and uq in the second line. The potential energy 
contribution is given by the particle density p times the 
dipolar energy per particle ~ P^P^- This gives in total 
the contribution ~ /i^p^ to the energy density. Expand- 
ing the density using p = po + up to second order in 
bp gives p"" ~ Po^p^ — Po(V(/)/7r)^. The clustered state 
of polar molecules has the same kinetic energy but a re- 
duced potential energy 

= 1^1 d:r[^X(Ve)2 + -|(V0n, 

as the potential energy is now given by the density of 
clusters pjq times the interaction energy of a single clus- 
ter [qpY /{q/ pY. This gives rise to the additional factor 
Comparing the Luttinger parameters in the g = 1 
and q > 1 cases gives 

1 Uq u 

uqKq = uK and ^ 1^ = t7 , 
q^ Kq K 

leading to 

K — qKo and u — u^/q . 

Therefore the Luttinger parameter K of the clustered 
Luttinger liquid is proportional to q and the decay expo- 
nent of the density-density correlation function is effec- 
tively suppressed by \/q, see discussion below Eq. ([2]). 
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